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Abstract 

The late-time growth of large scale structures is imprinted in the CMBR anisotropy through 
the Integrated Sachs Wolfe (ISW) effect. This is perceived to be a very important observational 
probe of dark energy. Future observations of redshifted 21-cm radiation from the cosmological 
neutral hydrogen (HI) distribution hold the potential of probing the large scale structure over a 
large redshift range. We have investigated the possibility of detecting the ISW through cross- 
correlations between the CMBR anisotropies and redshifted 21-cm observations. Assuming that 
the HI traces the dark matter, we find that the ISW-HI cross-correlation angular power spectrum 
at an angular multipole i is proportional to the dark matter power spectrum evaluated at the 
comoving wave number i/r, where r is the comoving distance to the redshift from which the HI 
signal originated. The amplitude of the cross-correlation signal depends on parameters related to 
the HI distribution and the growth of cosmological perturbations. However, the cross-correlation 
is extremely weak as compared to the CMBR anisotropies and the predicted HI signal. Even in an 
ideal situation, the cross-correlation signal is smaller than the cosmic variance and a statistically 
significant detection is not very likely. 
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I. INTRODUCTION 

In recent times, a host of independent observations, like Supernova-la galaxy 
surveys [3] and Cosmic Microwave Background (CMB) anisotropies 0, [sl, have indicated 
that the expansion of the Universe is accelerating Q|. This can be explained by a dark 
energy component, with an equation of state pj p = w{< — |). The cosmological constant, 
A, has emerged as a strong candidate for dark energy, as various observations 51] constrain 
w to be be close to —1. 

An indirect effect of A is that it causes a decay of the gravitational potential, when the 
universe evolves from the matter dominated to the the dark energy dominated era. This 
generates a weak anisotropy in the CMB temperature fluctuation, through the Integrated 
Sachs Wolfe (henceforth ISW) effect 71]. A non-flat spatial geometry would contribute to the 
ISW in the same way. However, CMB data largely constrain our universe to be spatially fiat 
(see ref. jsj) so, such effect of spatial curvature can be ignored in the first approximation. 
The late-time evolution of the gravitational potential is sensitive to the specific dark energy 
model. Therefore, the associated ISW anisotropy, can in principle be used to probe the 
nature of dark energy. 

It is difficult to separate the ISW signal from the primary CMB anisotropy, because it is 
intrinsically weak and it appears at large scales, where the error due to cosmic variance is 
large. 

Techniques to measure the ISW, use the cross-correlation of the CMB fluctuations, with 
fluctuations of some tracer of the large scale structure at a later redshift. Fluctuations in 
the primary CMB field and in the tracer are uncorrelated, so that this method allows one to 
single out the contribution solely due to the ISW. It is also important to note here, that the 
foregrounds and noise are not correlated between independent random fields. Recently ISW- 
large scale structure and ISW-weak lensing cross-correlations have been studied extensively 
(see ref. [isl]). These studies look at a median 2; ~ 1.5 and are in 

agreement with the ACDM model at 2cr ~ 3cr levels. 

Observations of redshifted 21 cm radiation of the spin-flip hyperfine transition from neu- 
tral hydrogen (HI) have the potential of probing the universe over a wide range of redshifts 
(200 > -2 > 0): from the dark ages to to the present epoch (eg. [16, 17, 18]). Recently, 



radio-optical cross-correlation study has detected a positive correlation between the opti- 
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cal gala xies (6dFGS) and HI fluctuations [19^. This vindicates the theoretical predictions 



21| about the possibility of using HI distribution statistically, as a probe of the large 



scale structure, without the need to resolve individual galaxies. Cross-correlation technique 
using the HI 2 1-cm radiation as one of the flelds has been considered for the study of cosmic 
reionization 



22 



23|,l2J,l25|,l26|. 



In this paper we study the use of diffused HI as a tracer of the large scale structure 
to probe dark energy induced ISW effect. We look at the cross correlation between the 
post-reionization (-2^6) fluctuations in the HI brightness temperature and the CMB. 

Redshifted 21 cm observations of neutral HI allow us to probe the universe as a function of 
redshift. The advantage of using HI tomography is that, we can probe the late-time cosmic 
history continuously over a range of redshifts. Radio telescopes (eg. currently functioning 
GMRT [48| and upcoming MWA [3] & LOFAR are aimed to map the large-scale 



distribution o 
details see 27 



HI at high redshifts. At redshifts < ^ < 3.5 we have f2gas ~ 10 ^ (for 



28 



291]). This implies that the mean neutral fraction of the hydrogen gas 



is xm = 50 Qgash'^{0.02/Qbh'^) = 2.45 x 10~^, which we assume is a constant over the the 
entire redshift range z < 6. 

The redshifted 21 cm radiation seen in emission in this redshift range, from individual 
clouds is rather weak (< lOyuJy). This makes its detectability dubious, with existing ob- 
servational facilities. (There might be considerable magniflcation caused by gravitational 
lensing which may enhance detection chances). Statistical distribution of HI however 
produces a weak background in radio observations. This radiation has the information about 



the HI fluctuations in probed redshift range 20|, l2ll]. CMB map of a large portion of the 



sky and a corresponding HI map would allow us to compute the cross-correlation power 
spectrum and hence independently quantify the cosmic history at redshifts z < 6. 

II. FORMULATION 



The CMB brightness temperature fluctuation along the direction of the unit vector n is 



described by 



31 



32| 



AT(n) = T <{ ( + V ■ n + $ ) + H dr][^ + ^] 

^ / LSS JrjhSS 



(1) 
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where T is the CMB temperature at present. Here, under the assumption of instantaneous 
recombination, the Sachs Wolfe effect (first term) is evaluated at the last scattering surface 
(LSS) and the ISW effect (second term) is integrated from the LSS to the present epoch. 
The scalar potentials $ and ^ are the metric perturbations in the conformal Newtonian 



gauge 
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3J], the dots refer to differentiation with respect to the conformal time r] and we 



shall use r = rjo — t] to denote the comoving distance to the conformal time rj. 



In the absence of anisotropic stress we have $ = 



33| and the ISW term is 



rvo 

Ar(n)^^^ = 2r/ dri^{rh,ri). (2) 
Expanding this in the basis of spherical harmonics 

oo 

AT{nY'^ = J2af::fYUn) (3) 

£,m 

and using the identity 

47rY,HYMkr)Y;^{mUn) (4) 



i.m 



we have 



a 



ISW 



/ d^^k,^)jeikr)Y;jk) (5) 



where $(k, r/) is the Fourier transform of $(r, ■?]), and je{x) is the spherical Bessel function. 

For sufficiently sub-horizon scales the gravitational potential can be related to the matter 
density fluctuations 6 via the Poisson equation. In Fourier space this takes the form 

Further, retaining only the growing mode of density perturbations 6(k,a) = 5(k)D+(a) 
we have 

i,(k,^) = i^^^|.(k,r/) (7) 
a 

where 

d{lna) 

which we use in eq. ([5]) to calculate ■ 
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The HI 21-cm brightness temperature fluctuations from redshift zhi can, in Fourier space, 



be written as 



35| 

AH/(k)=fxH/(& + //i')5(k,a) (9) 



where xhi is the mean HI fraction, /i = k ■ fi and 



f(.).,0n,K(l..,^(5^)(Y)^ (10) 

Here it has been assumed that the HI traces the underlying dark matter distribution with a 
possible bias h. On the large scales under consideration, where the matter fluctuations are in 



the linear regime, it is reasonable to assume that the baryonic matter follows the under 



ymg 
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dark matter distribution. The term //i^ has its origin in the HI peculiar velocities 
which have also been assumed to be caused by the dark matter fluctuations. It should be 
noted that all the terms on the rhs. of equation (l9l) are to be evaluated at the redshift zhi at 
which the HI signal originated. Note that one should include a normalized window function 
W{z) in eq. ([9]) describing the spectral response of an instrument |25i]. On scales of our 
interest {^ < 100), the spectral resolution of the instrument can however be assumed to be 
much smaller than the features in the HI signal and W(z) can be approximated by a 
Dirac delta function, so that eq. (I9l) is, a reasonably good approximation. 

Expanding the HI signal in terms of spherical harmonics and proceeding as before we get 

^<5(k,a)T,(fcr)F4(k) (11) 

where 

I,{x)=h3,{x) - f'^ . (12) 

We use equations ([5]) and ( ITT]) to calculate Cf^~'^^ the cross correlation angular power 
spectrum between the HI 21 cm brightness temperature signal and the CMBR ISW signal 
deflned through 

Note that Cf^~^^^ also depends on z^j the redshift from which the HI signal originates, or 
equivalently on i> = 1420MHz/(l + zhi) the frequency of the HI observations, but we do 
not show this explicitly here. We obtain 



/rvo 
dk P{k)UkrHi) / dvF{7^)je{kr 



(14) 



where P{k) is the present day dark matter power spectrum, 



Aiz) = -fiz)xHiD4z)^^^^ (15) 



and 



FM ^ £i(i^msi (16) 

For large i we can use the Limber approximation [9|, 137[ which allows us to replace the 
spherical Bessel functions by a Dirac deltas 6d{x) 



Mkr)^^l^^6n{i + l-kr) (17) 
whereby the angular cross-correlation power spectrum takes the simple form 

« ^-^^P(^) (18) 

where P{k) is the present day dark matter power spectrum and all the other terms on the 
rhs. are evaluated at zhi- 

we also have. fo. co,„pa.so.. .he HI-HI a^gula. power speet™,„ (A.) 0. which 
describes the statistical properties of HI fluctuations at two redshifts (corresponding fre- 
quencies being u and u + Au). Using the 'flat sky' approximation 38|] we have C^^^Au) is 
given by 

rp2 roc 

C^\Au) = — .xli/^i / dkn cos(A;|| < Au) \b + ff?] ' P{k) (19) 
T^^i Jo 

where r is the comoving distance corresponding to the redshift zhi or equivalently frequency 



z/ = 1420MHz/(l + zm), K = drjdu and k = ^kj + (//r)^. 

The function C^^Au) is a direct observational estimator of the HI fluctuations at redshift 
Zm- This does not require us to assume an underlying cosmo logical model (eg. [391]). Here we 
have taken the special case where Au = and we shall henceforth refer to the corresponding 
power spectrum as C^^ . 



III. RESULTS 

Figure [1] and Figure [2] respectively show the predicted HI-ISW cross-correlation an- 
gular power spectrum Cf^~^^^ and the HI-HI angular power spectrum C^^ for a few 
redshifts zhi in the range 0.5 < zhi < 5. We have used equations [14] and [19] to 
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FIG. 1: The HI-ISW angular power spectrum for redshifts z = 0.5.1.0,2.0,3.0,4.0,5.0 (top to 
bottom). 
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FIG. 2: The HI angular power spectrum C^' at redshifts z = 0.5.1.0,2.0,3.0,4.0,5.0 (top to 
bottom). 



calculate the cross-correlation angular power spectrum and HI power spectrum respec- 
tively. The approximated equation f|T8l) is useful for qualitative description of the re- 
sults. We have assumed the currently favored ACDM cosmological model with parameters 
(n„o,^Ao,/i,(^8,^s) = (0.28,0.72,0.7,0.82,0.97) \^,^. The bias, b for the post reionization 
HI on large scales is assumed to be linear. We have taken 6 = 1 as the fiducial model. 
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However, it is important to note that HI in the post reionization epoch is assumed to be 
distributed in high column density clouds which could be more biased with respect to the 
underlying cold dark matter distribution. 

The shape {i dependence) of the cross-correlation signal reflects the shape of the matter 
power spectrum P{k) (eq. [TS]). We flnd a peak in C^^^^^^ at i = rkeq, where keq is the 
wave vector corresponding to the matter radiation equality. For different redshifts zhi the i 
value corresponding to this peak scales as £ oc r, the comoving distance to the redshift zhi- 

The amplitude of the cross-correlation signal Cf^~^^^ depends on a product of various 
terms some of which (T, xhi-, b) depend on the HI distribution and others (-D+, /, H) which 
depend on the cosmological model. The dimensionless term / quantifles the growth of the 
dark matter perturbations, and the ISW effect is proportional to / — 1. We have / = 1 in 
cosmological models with no dark energy, and we do not expect to have any ISW effect in 
such models. The term / — 1 is a sensitive probe of dark energy. The amplitude of Cf^'^^'^ 
contains this information combined with unknown parameters related to the HI distribution. 
It has been recently proposed that observations of the HI fluctuations at low z can be used to 
estimate cosmological parameters 4o|, 41|- It is in principle possible to combine observations 
of Cf^~^^^ and Cf^ to jointly estimate parameters of the HI distribution and the background 
cosmological model. 

IV. DETECTABILITY AND CONCLUSIONS 

Here we estimate the viability of detecting the HI-ISW cross-correlation signal. The 
cosmological HI signal is weak and buried under the foregrounds which are orders of mag- 



nitude higher than the signal 



16 



38 



39 



42 



43| . This is a serious observational problem 



for auto-correlation studies involving the 21cm radiation. One may separate the foreground 
components by noting that HI signal (a line emission) decorrelates beyond a certain fre- 
quency separation whereas the foregrounds remain correlated over large frequency separa- 
tions. We shall subsequently assume that foregrounds have been removed. Moreover, the 
cross-correlation signal is less affected by foregrounds and other systematics. This is be- 
cause, many of the foregrounds and noise are expected to be uncorrelated between the two 
maps. 

The uncertainty in estimating the cross correlation signal is the sum, in quadrature, of 
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FIG. 3: The cumulative S/N by collapsing all multiploes less than £ for different redshifts. 



the instrumental noise and the cosmic variance. While the system noise can, in principle, be 
reduced by increasing the duration of the observation the cosmic variance sets a fundamental 
limit in deciding whether the signal can at all be detected or not. 

The cosmic variance of the cross-correlation angular power spectrum Cf^~^^^ is {q] 

riCMBr^Hi 



a 



(20) 



(2£+l)v^/,A£ 

where Cf^^^ is the CMB angular spectrum for which we have used the WMAP5 results 
5l| . A£ is the width of bands in £ and fs is the fraction of the sky common to both 
the CMBR and HI observations. We have used M = 10 for i < 100 and Ai = 100 for 
i > 100, and have assumed the most optimistic possibility fs = I for our estimates at 
redshift ^hi = 0.5. Different frequencies channels in the frequency band of HI observations 
provide A^^^ independent estimates of the HI signal which cause a reduction in the cosmic 
variance by a factor l/^/Wl- Here we have assumed that the HI observations are carried 
out across a bandwidth of 32 MHz centered around zm = 0.5 and the HI signal is assumed 
to be independent at frequency separations of ~ IMHz |21|], which gives Nc = 32. Using 
these to estimate the signal to noise ratio S/N = Cj^^~^^^ /a we find that S/N < 0.45 for all 
zhi and i and a statistically significant detection is not possible in such cases. It is possible 
to increase S/N collapsing the signal at different multipoles i. To test if a statistically 
significant detection is thus feasible we have collapsed all multipoles less than i to evaluate 
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FIG. 4: The cumulative S/N on combining data upto a certain redshift. 



the cumulative S/N defined as [221. |44| 



S 



(2£+l)v^/,(C, 



HI-ISW\ 2 



(jCMBfjHI 



(21) 



Results are shown in Figure [3] for various redshifts (0.4 < z < 3). We find that the contri- 
bution in the cumulated S/N comes from £ < 50 at all redshifts that we have considered. 
The cross-correlation signal is largest at {z ~ 0.4) and is negligible for {z > ?>). We further 
find that although there is an increase in S/N on collapsing the multipoles it is still less 
than unity. This implies that a statistically significant detection is still not possible. Thus, 
probing a thin shell of HI doesn't allow us to detect a cross correlation, the signal being 
limited by the cosmic variance. 

21 cm observations have the advantage that one may probe various redshifts by tuning 
the frequency of radio observations. This enables us to optimally combine the signal from 
a large number of thin shells over a continuous range of redshifts. We have considered a 
range of redshifts (0.4 < ;z < 3 or 1000 > v > 350) and combined the signal for independent 
observations at 32 MHz separations in this range. The S/N cumulated upto a certain 
redshift is shown in Figure HI This indicates an increase in the S/N. A cumulated S/N 
of ~ 1.6 is attained for redshift upto z = 2 and there is hardly any increase in S/N on 
cumulating beyond this redshift. This is reasonable because the contribution from the ISW 
effect becomes smaller beyond the redshift z > 2. This S/N is the theoretically calculated 
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value for an ideal situation and is unattainable for most practical purposes. Incomplete 
sky coverage, and foreground removal issues would actually reduce the S/N and attaining 
a statistically significant level is not feasible. We conclude that, within the paradigm of 
ACDM cosmology, though there is a weak positive correlation between the CMBR ISW and 



HI, the signal is much weaker than the individua 
unlikely. Certain modified gravity models (eg. 45 



auto-correlations and a detection is quite 



46|) may allow the quantities (-D+, /, H) 



to be different from what they are in the ACDM model (considered here) 47| and may lead 
to an increase of the S/N. However, since the cross-correlation signal is significant only at 
large scales we don't expect the S/N to be much different from the ACDM predictions. 
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